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Abstract 

We discuss the appearance of spurious solutions of few-body equations for 
Faddeev amplitudes. The identification of spurious states, i.e., states that lack 
the symmetry required for solutions of the Schrodinger equation, as well as 
the symmetrization of the Faddeev equations is investigated. As an example, 
systems of three and four electrons, bound in a harmonic-oscillator potential 
and interacting by the Coulomb potential, are presented. 



I. INTRODUCTION 

One of the most viable approaches for solving the few-body problem is the Faddeev 
method It has been successfully applied to solve the three- nucleon bound-state problem 
for various nucleon-nucleon potentials The most complex calculations of this kind 

include up to 82 channels, when all the j < 6 waves are taken into account [|] . 

The Hamiltonian in the Schrodinger equation is Hermitian and the solutions for a system 
of fermions, for example, are antisymmmetrized. On the other hand, the Faddeev equations 
are non- Hermitian and the trial wave functions used for a system of fermions, e.g., are not 
fully antisymmetrized. Also, there are three Faddeev equations for the three-body system 
acting on the same variables as in the Schrodinger equation. One, therefore, could expect 
that there should be three times as many solutions as for the Schrodinger equation. Then, it 
is not surprising that spurious solutions of the Faddeev equations exist. A spurious solution 
is an eigenstate, which does not have the symmetry required by the Schrodinger equation, 
e.g., antisymmetry for a system of identical fermions. A spurious component in a solution 
of the Faddeev equations was found analytically for the first time for the ground-state of 
three identical particles bound in the harmonic-oscillator (HO) potential by J. L. Friar et 
al. 0. Spurious components of the Faddeev amplitudes were then observed also for the 
excited states in extensions of this work for three identical particles M, three nonidentical 
particles ||, as well as for four identical spinless particles |L0| . Spurious solutions of Faddeev- 
like equations were investigated by Evans and Hoffman |IT|, and the existence of spurious 
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solutions of Faddeev equations for three identical particles was recently demonstrated by 
Rudnev and Yakovlev 



12 . In Ref. [12 



such solutions were constructed for states of zero 
total angular momentum. In addition, spurious solutions of the Faddeev equations for three 
nonidentical particles interacting by central potentials were investigated by V. V. Pupyshev 
fl~3|j . In Ref. |TJ|], we noted the appearance of spurious solutions, whose number exceeded 
the number of physical solutions, of the Faddeev equations for a three-nucleon system solved 
in an HO basis in a shell-model approach. Also, a spurious solution was reported recently 
in a three-body model calculation of ^Be |15| . 

In the present paper we investigate the appearance of spurious states in the Faddeev 
formalism and their identification in a systematic manner. We use three- and four-electron 
systems bound in an HO potential as an example for illustration and quantification of the 
problem. 

In section [H] we discuss the three fermion system. A generalization for the four-body 
system is presented in section |T|. Conclusions are given in section |V[ 



II. THREE-BODY SYSTEM 

Our discussion is quite general. However, we prefer to illustrate our points by using a 
particular simple example, namely, a system of electrons bound by an HO potential and 
interacting by the Coulomb potential. We consider, therefore, the following Hamiltonian 



.4 



i=l I 



r? 1 
2m 2 



A 



EMri-rj). (1) 



Eigenstates of the Hamiltonian ([[]) are antisymmetric with respect to the exchange of any 
electron pair. 

In the Faddeev formalism for a system of three identical particles, i.e., A = 3, the 
following transformation of the coordinates 

Un-r 2 ), (2a) 



2rl 

3 L2 



[Un+f 2 )-r 3 ] , (2b) 



and, similarly, of the momenta, is introduced. It brings the one-body HO Hamiltonian from 
Eq. (HD into the form 



Ho = ^ + + jf- + Imtff , (3) 

2m 2 2m 2 

with the trivial center-of-mass term omitted. Eigenstates of H , 

\nlsj,Af£J,J) , (4) 

can be used as the basis for the Faddeev calculation. Here n, I and A/", C are the HO 
quantum numbers corresponding to the harmonic oscillators associated with the coordinates 
and momenta r, p and y, q, respectively. The quantum numbers s, j describe the spin and 
angular momentum of the relative-coordinate partial channel of particles 1 and 2, J is the 
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angular momentum of the third particle relative to the center of mass of particles 1 and 
2 and J is the total angular momentum. It is obvious that the basis @) (similarly for 
trial states used in traditional Faddeev calculations), while antisymmetrized with respect 
the exchange of particles 1 «-> 2, is not antisymmetrized with respect the exchanges of 
particles 1 •*->• 3 and 2^3. On the other hand, the physical solutions corresponding to the 
solutions of the Schrodinger equation must be totally antisymmetrized. For our example 
system, the states (^) form a complete orthonormal basis. The physical solutions are such 
linear combinations of the basis states that have the proper antisymmetry for all the particle 
exchanges. On the other hand, one must expect that there are more basis states (Q) than the 
total possible number of linearly independent antisymmetrized linear combinations of the 
states (Hj). Consequently, it is natural that spurious states will appear, when the calculation 
is performed using a basis of the type (f|). 

The first Faddeev equation can be written in the differential form as 

{H + V 3 - £)0 3 = - Wi + <h) , (5) 

where V 3 is the two-body interaction between particles 1 and 2 and 4>i,4> 2 , and 3 are the 
three Faddeev components. The other two Faddeev equations are cyclical versions of the 
Eq. (H). Equation (|5]) can be re- written in the form 

H\1>) K = E\1>) K , (6) 

with 



H = H + V{f)T . 



(7) 



Here, V(r) is the interaction between particles 1 and 2, e.g., V3 of Eq. or Vc(y2r) of 
Eq. (U); \i/j)k corresponds to the Faddeev amplitude 3 in Eq. ([5]); and T, which has the 
properties of a metric operator P,16|, is given by 



r = 1 + r ( - } + r (+) 



(8) 



with 7"( + ) and T* - ) the cyclic and the anticyclic permutation operators, respectively. 

For the basis (|4]), which we use as an example in the present paper, we obtain a formula 
for the matrix elements of T* - ) + by simplification of the expression (10) in Ref. []T4 
namely 



(niZiSaji,MAJi, J\T { - ] + T {+) \n 2 l2S232^ 2 C 2 J 2 , J) = 5 Nl 

h si jx ) ( h s 2 j 2 



LS 



L S J 



L S J 



A ? 2 

1 I « 

2 2 61 
\ S S 2 



S1+S2— C%— I 



(9) 



where N t = 2n^ + k + 2A^ + C t ,i = 1,2, ] = s/Tj + 1 and (M 1 C 1 n 1 l 1 L\n 2 l 2 J\f 2 C 2 L) 3 is 
the general HO bracket for two particles with mass ratio 3 as defined, e.g., in Ref. ||17|| . 
Similar expressions for different bases are described, e.g., in Refs. [TROll. It follows from 
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the symmetry properties of the states (|4]) and of the HO brackets that the contributions of 
and in @ are identical. 
We note that the eigensystem of the operator T (||) consists of two subspaces. The first 
subspace has eigenstates with the eigenvalue 3, which form totally antisymmetric physical 
states, while the second subspace has eigenstates with the eigenvalue 0, which form a not 
completely antisymmetric, unphysical subspace of states. Although we found these proper- 
ties of T by direct calculation using the relation (01), it is, in fact, a general result. The same 



structure of eigenstates was also reported in Ref. |12| using a different basis. The eigenvalue 
structure follows from the fact that |T has the properies of a projection operator. 

The Hamiltonian H (|7]) is non- Hermit ian. By solving the equation (|6]) one obtains the 
right (ket) eigenstates, while by acting with H to the left, one gets the bi-orthogonal left (bra) 
eigenstates. In the basis with physical and spurious states separated, e.g., those obtained 
from the diagonalization of T, the Hamiltonian matrix takes the form 

/(Ph|#|Ph) \ 

^ (s P |#|Ph) (Sp|#|s P > ; ■ {U> 

Formally we obtain the Hamiltonian matrix (|It]) from the Hamiltonian matrix in the ba- 
sis of the type (||) by an orthogonal transformation that transforms the basis (|]) into the 
eigenstates of the operator T (RJ). In fact, this situation is very much analogous to the 



well-known properties of the Dyson boson- mapped systems pQj . It follows from the struc- 
ture of the Hamiltonian matrix fliCf ) that its right spurious eigenstates Sp)k are not 
contaminated by the physical states. Similarly, the left physical eigenstates b("0; Ph| are 
not contaminated by the spurious ones. On the other hand, its right physical eigenstates 
|-0; Ph)x may have spurious admixtures; similarly the left spurious eigenstates b(0; Sp| may 
have physical admixtures. As a right spurious eigenstate has no physical admixtures, it must 
be annihilated by the action of the operator T (|8|), e.g., 

T|^Sp> K = 0. (11) 

The action of the opeartor T on the right eigenstates, thus, serves as a test for identification 
of spurious states among the calculated eigenstates, in analogy to identification of spurious 
states in the Dyson boson mapping by the means of the (so-called) 1Z operator The 
relation ([□]) can be expressed in terms of the Faddeev amplitudes appearing in Eq. ([|) 
as 0i + 02 + 03 = 0. We note that the sum of the three Faddeev components has to 
be indentified with the solution of the Schrodinger equation. If one finds solutions of the 
Faddeev equations, which are unphysical, e.g., that do not have the proper symmetry with 
respect to the particle exchanges, the corresponding three Faddeev components have to add 
up to zero. This is the only possibilty to avoid a contradiction. 

It is possible to avoid the spurious state problem completely by hermitizing the Hamil- 
tonian (0) on the physical subspace, where it is quasi-Hermitian (see the discussion of quasi- 



Hermitian operators, e.g., in Ref. |L6]]). The Hermitized Hamiltonian takes the form 

H = H + T l / 2 V(r)T 1 / 2 , (12) 
where T operates on the physical subspace only, e.g., it has only the eigenvalues 3. 
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Let us mention several relations between the eigenstates of the Faddeev Hamiltonian 
H (ffj) and the symmetrized Hamiltonian H (12). First, we have for the left and right 
eigenstates of H 



B <V;Ph|= K<^;Ph|T. (13) 

The eigenvectors \$) of the Hamiltonian H flT2|) are related to the eigenvectors of the Hamil- 
tonian H (0) by 

|M/) = T 1 /V;Ph)K, (14a) 
(^1 = b(^; Ph|T~ 1/2 . (14b) 

Note that we write explicitly T instead of T only when an inversion is needed. For matrix 
elements of a general operator O we then have 

(* f \o\%) = B (^/; Phir-^or 1 / 2 !^; s P ) K . (15) 

In addition to the relation (|TT|), there is another possibility of identification of spurious 
states. It follows from the properties of the matrix elements of a general operator O between 
the physical and spurious eigenstates of if (|7|), namely 

B (V/; PhlT-^OT 1 / 2 ^; Sp) K = , (16a) 

B <V/; Spir-^or 1 / 2 !^; Ph) K ^ o . (16b) 

These relations follow from the properties of the left and right eigenstates as discussed in 
the paragraph after Eq. flTOl). For an operator that commutes with T the square roots of 
T may be omitted in Eqs. QT5| , |TB[ ) . Also, in the above relations the substitution J2i<j Oy = 
T 1 ^ 2 0{f)T 1 / 2 can be used for a two-body operator Oij depending on (fj — fj). 

Let us return to our specific system example. The metric T (j^) is diagonal in N = 
2n + 1 + 2JV+ C. It follows from the expression (|9]) that any basis truncation other than one 
of the type iV < iV max would lead, in general, to mixing of physical and unphysical states. 
Here, iV max characterizes the maximum of total allowed harmonic-oscillator quanta in the 
basis. At the same time, the truncation into total allowed oscillator quanta iV < iV max 
preserves the equivalence of the Hamiltonians (0) and ( |T2"D on the physical subspace. A 
general consequence of this observation is that an improper truncation in the treatment of 
T leads to physical and spurious state mixing. 

In Table | we present the dimensions (D) of the basis (|J) corresponding to a particular 
iV = 2n+l+2Af+£ together with the number of physical states (Ph) and spurious states (Sp) 
for the three-electron system with J 71 = l/2~. Apparently, the number of physical states is 
about a third of all basis states. We also present the ground state and the first excited state 
energies obtained with the basis restricted by iV < iV max , where iV max corresponds to the 
number in the first row. We used the HO energy HQ = 0.5 atomic units (a.u.). The physical 
eigenenergies are shown without the trivial center-of-mass contribution. It is immediately 
seen from Eqs. ([?]) and (PD that the spurious states have the energies corresponding to the 
unperturbed Hamiltonian H (|3p. Typically in the search for a bound state in the Faddeev 
calculations the lowest state would be physical and the spurious state existence would be 
unnoticed. The present electron system with a repulsive Coulomb interaction is interesting 
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because of the fact that the physical ground state is the tenth state as can be deduced from 
the Table |. There are 2 spurious states with the unperturbed HO energy 2 a.u. and 7 
spurious states with the energy 3 a.u. We note that for the discussed system the ground- 



state energy obtained by the Stochastic Variational Method (SVM) [[22|] is 3.26324 a.u. |23 
(after subtracting 0.75 a.u. for the center-of-mass energy). In our HO basis calculation two- 
decimal place precision is obtained rather rapidly. A further improvement of the precision 
is, however, slow. We performed calculations up to N max = 39, where we obtained the 
ground state energy of 3.2634 a.u. and 4.0446 a.u. for the first-excited state. A substantial 
acceleration of convergency can be achieved by employing the effective interaction approach 
in a manner similar to that discussed in Ref. [Q. By replacing the interaction V(r) in Eqs. 
(0.|3) by V e ff(f) we reach the SVM ground-state result for iV max = 27 and the for the first 
excited state we then obtain 4.04458 a.u. 



III. FOUR-BODY SYSTEM 

To demonstrate that the problems discussed prevail also for the Faddeev-type approach 
to systems with more than three particles, we present briefly the extension of the studied 
system to four electrons. We use the Hamiltonian ([!]) with A = 4. By introducing the 
coordinate (and momenta) transformation 



r = JUn-r 2 ), (17a) 



y=^IWi+r 2 )-?3]i (17b) 
z = &[\(r 1 +r 2 + r 3 )-rt], (17c) 

the one-body part of the Hamiltonian ([]]) is obtained as 

H = ^ + -m^V + ^- + -mQ 2 f + — + -mtfz* , (18) 
u 2m 2 2m 2 y 2m 2 ' v ; 

with the trivial center-of-mass term omitted. 

A possible generalization of the Faddeev equation (§) for four identical particles can be 
written in the form 

W(123)4> = W(123)4> , (19) 

with 

##(123)4) = #(#(123)4) + \(Vi2 + V 13 + V23X |^(123)4) + l'0(432)l) + 1^(134)2) + 1^(142)3)) , 

(20) 

and 

(1^(123)4) + |V>(432)l) + 1^(134)2) + 1^(142)3)) = (1 - T U - T 24 - 2^)1^1(123)4) = 24^(123)4) • 

(21) 
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Here, | '0(123)4 ) is a four-fermion Faddeev amplitude completely antisymmetrized for particles 
1,2, and 3. There are three other equations that can be obtained from Eq. (|i~9l) by permut- 
ing particle 4 with particles 1, 2, and 3. Their sum then leads to the Schrodinger equation. 
We note that the present equations are different from the traditional Faddeev-Yakubovsky 



equations which combine Faddeev amplitudes depending on two sets of Jacobi coordi- 
nates. As we are working with a complete orthonormal basis, it is sufficient and convenient 
to use a single set of Jacobi coordinates defined by the relations (O). Unlike the Faddeev 
amplitudes used typically in the Faddeev-Yakubovsky equations, the amplitudes appearing 
in Eq. (|TjJ) are antisymmetrized with respect to the first three particles. The present equa- 
tions allow to employ easily three-body interactions or three-body effective interactions. The 
latter property makes them particularly useful for an extension of shell-model calculations 
for four nucleons p5[ . 

The starting point for the present four-electron calculation is the basis 



\NxiJx,n g l x Ji, J) 



(22) 



with the three-fermion part given by the antisymmetrized eigenstates of T (^) corresponding 
to eigenvalue 3, e.g., 

WiUi) = £ <8&£*\rd8j,M£J*, Ji) , (23) 

where Aq = 2n + I + 2Af + C and i counts the eigenstates of T with the eigenvalue 3 for 
given Ni and 3\. 

As for the three-particle transposition operators @, a compact formula can be derived 
for the matrix elements of the four-particle transposition operators in the basis (0), e.g., 

(Ail«l JiL, Tlzhhh^h, J|7l4 + ^24 + AirZr Jir, TI z rI z rJ^ J) 
_ _X \^ „N 1L i L J lh „N 1R i R J 1R f2 f2 r.2 q2 f 2 q2 

- °N h ,N R 2^ C n L / LSL j L Ar L £ L J3 L C n R ZR SR j R Ar R £ H ,J3 R ^lL^lR D lL D lR^2 D 2 

\ SR S 



XjLkj 3 LJ3Rj4LJARJlLJlR{-l) SlL+Sm * 2 

Llh S'lL JlL 



1R 

2 $2 S'lL 



X 



^R S R Jr. 
Al o «^3R 



XL 



12 



k s L j L 

Llh Sil JlL J [ JlR Sir JlR 

l' \ Ir L2 L' 1 J Zr L 2 J' 
I ^R £r JlR J \ ^L ^ JlL 

'zR+^lL f /" 1 \lzL 



hh \ <Jih 
L2 S2 J 



L\r Sir Jir 

IzR 2 ^ 4R 
J2 S2 J 



-I; 



ShSR< f 



1 1 c 

2 2 SR 



SlL Sl 



X (_l)'«a+ML ^(-l)^ L (n7 , n 2L / zL L> Z RUAr R £RL08^L/LAr L /: L L 1L |n7 , nR/RL 1L )3 
+(-l)'a-*+ £ a-fc-^<n* L Z^ 



+^,i K ^s h ,s I{ S A f h:nl 5c h ,i'(-l 



(^ Z L^LA' L /: L L / |A'R£Rn 2 R/ 2 RL / ) f 



(24) 



where A x = 2n x + Zx + 2A/x + £x + 2n 2X + l z x = N 1X + 2n zX + l z x, X = L, R. We will 
give more details on the derivation, symmetry properties, as well as a generalization that 
includes the isospin quantum numbers, elsewhere [2~E . 

As the structure of the equation ( p!9|) is the same as the Faddeev equation (|6]), the 
discussion of the spurious state problem for the three-fermion system can be extended to 
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the four-fermion system as well. In particular, eigenstates of the operator defined by the 
relation ( plj) can be subdivided into two subspaces. A physical subspace corresponding to 
the eigenvalue 4 spanned by totally antisymmetric states and a spurious subspace spanned 
by eigenvectors corresponding to the eigenvalue 0. In the basis of the eigenstates, the 
Hamiltonian H appearing in Eq. ([[£]) has the same structure as shown in Eq. (|T0|). Its 
spurious eigenstates can be identified by the action of T± on the right eigenstates, like in Eq. 
(0), e.g., T 4 |^ (123 )4; Sp) K = 0. 

It is possible to symmetrize the Hamiltonian H on the physical subspace and eliminate 
the spurious state problem at the same time. The symmetrized Hamiltonian then takes the 
form 

H = H + % 1/2 \(V 12 + V 13 + V 23 )% 1/2 , (25) 

where operates only on the physical subspace. 

In Table |T|we show the dimensions (D) of the basis (p2"|) corresponding to particular N = 
Ni + 2n z +l z together with the number of physical states (Ph) and spurious states (Sp) for the 
four-electron system with J 71 = + . The relative number of physical states decreases with N 
and approaches about a fourth of all basis states for larger N. In the last two rows of Table [D 
we present the ground state and the first excited state energies obtained in the calculations 
with the basis restricted by iV < iV max , where iV max corresponds to the number in the first 
row. As in the three-electron calculations the physical eigenenergies, which were obtained 
with the HO energy hQ = 0.5 a.u., are shown without the trivial center-of-mass contribution. 
We note that the studied system, when described by Eq. (|T9|) , has a spurious ground state 
with the energy 3.25 a.u., corresponding to the unperturbed HO Hamiltonian H (ITH). We 



note that the preliminary result of the ground-state energy of the studied system obtained 
by the SVM is 5.6002 [|3|] (after subtracting 0.75 a.u. for the center-of-mass energy). As for 
the three-electron calculations, we get a fast convergence to two-decimal place precision in 
our HO basis calculation. A further improvement of the precision is slow and requires larger 
model spaces. As in the three-electron calculation, a substantial acceleration of convergence 
can be achieved by employing the effective interaction approach [[y]]. By replacing the 



interaction V by a two-body effective interaction, we obtained a ground-state energy of 
5.5991 a.u. and a first-excited-state energy of 5.6958 a.u. for iV max = 14. An additional 
speeding up of the rate of convergence can be achieved by using the three-body effective 



interaction 25 



IV. CONCLUSIONS 

To summarize, we have investigated the appearance of spurious states in the Faddeev 
formalism for few-body systems. Depending on the studied system, such states may appear 
even among the lowest eigenstates. We have also discussed how the spurious states can be 
identified and how they can be eliminated by symmetrization of the Faddeev equations. We 
also noted that any improper truncation in the treatment of the particle exchange operators 
(e.g., T, I4) may lead, in general, to the mixing of physical and spurious states. 

We used three- and four-electron systems bound in an HO potential as an example. Due 
to the repulsive character of the Coulomb interaction, such systems have spurious ground 



S 



states, when solved in the Faddeev formalism. By examining the basis dimensions, we have 
illustrated that the number of physical states for the three-electron system is about a third 
of the total number of basis states. (It is exactly one-third, when the total number of basis 
states is divisible by three). For the four-electron system we formulated equations for the 
Faddeev amplitudes antisymmetrized for the first three particles. Using this approach we 
observed that almost 3/4 of all the states were spurious. 
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TABLES 



N 


1 


3 


5 


7 


9 


11 


13 


15 


17 


19 


21 


D 


3 


10 


21 


36 


55 


78 


105 


136 


171 


210 


253 


Sp 


2 


7 


14 


24 


37 


52 


70 


91 


114 


140 


169 


Ph 


1 


3 


7 


12 


18 


26 


35 


45 


57 


70 


84 


Ei 
E2 


3.4105 


3.2787 
4.1002 


3.2699 
4.0497 


3.2669 
4.0472 


3.2656 
4.0459 


3.2649 
4.0453 


3.2645 
4.0450 


3.2643 
4.0449 


3.2641 
4.0448 


3.2639 
4.0447 


3.2638 
4.0447 



TABLE I. Results for three-electron system bound in a HO potential are presented. The first 
line shows the total of HO quanta corresponding to the relative motion in the basis states. The 
dimensions corresponding to the state J n = l/2~ are shown in lines 2-4. The line 2 displays the 
full dimension (D), while the lines 3 and 4 present the number of spurious (Sp) and physical (Ph) 
states, respectively. The lines 5 and 6 show calculated ground-state and first-excited-state energies, 
respectively. The calculations were done for h£l = 0.5 a.u. The ground-state center-of-mass energy 
(0.75 a.u.) is not included in the energies shown. Energies corresponding to a particular N were 
obtained in calculations, where all physical basis states up to N were included. 



N 





2 


4 


6 


8 


10 


12 


14 


D 





4 


21 


68 


178 


391 


767 


1390 


Sp 





2 


14 


48 


128 


284 


564 


1024 


Ph 





2 


7 


20 


50 


107 


203 


366 


Ei 




5.8829 


5.6346 


5.6094 


5.6043 


5.6022 


5.6012 


5.6006 


E2 




6.0239 


5.7396 


5.7116 


5.7034 


5.7003 


5.6988 


5.6979 



TABLE II. Results for four-electron system bound in a HO potential are presented. The first 
line shows the total of HO quanta corresponding to the relative motion in the basis states. The 
dimensions corresponding to the state J w = + are shown in lines 2-4. The line 2 displays the 
full dimension (D), while the lines 3 and 4 present the number of spurious (Sp) and physical 
(Ph) states, respectively. The lines 5 and 6 shows calculated ground-state and first-excited-state 
energies, respectively. The calculations were done for Ti£l = 0.5 atomic units. The ground-state 
center-of-mass energy (0.75 a.u.) is not included in the energies shown. Energies corresponding to 
a particular N were obtained in calculations where all physical basis states up to N were included. 
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